Difficulties for Compact Composite Object Dark Matter 
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It has been suggested "that DM particles are strongly interacting composite macroscopically 
large objects ... made of well known light quarks (or ... antiquarks)." In doing so it is argued 
that these compact composite objects (CCOs) are "natural explanations of many observed data, 
such as [the] 511 keV line from the bulge of our galaxy" observed by INTEGRAL and the excess 
of diffuse gamma-rays in the 1-20 MeV band observed by COMPTEL 17\. Here we argue that the 
atmospheres of positrons that surround CCOs composed of di-antiquark pairs in the favoured Colour- 
Flavour-Locked superconducting state are sufficiently dense as to stringently limit the penetration of 
interstellar electrons incident upon them, resulting in an extreme suppression of previously estimated 
rates of positronium, and hence the flux of 511 keV photons resulting from their decays, and also in 
the rate of direct electron-positron annihilations, which yield the MeV photons proposed to explain 
the 1-20 MeV excess. We also demonstrate that even if a fraction of positrons somehow penetrated 
to the surface of the CCOs, the extremely strong electric fields generated from the bulk antiquark 
matter would result in the destruction of positronium atoms long before they decay. 

PACS numbers: 95.35.-(-d 



The conventional view is that dark matter is both cold 
- non-relativistic - and so weakly interacting that it is 
collisionless. Its dynamical evolution within astrophysi- 
cal systems is then governed entirely by gravity. Canon- 
ical particle physics models of dark matter yield inter- 
actions with ordinary and dark matter that are indeed 
irrelevant at current astrophsyical densities and energies. 
This conventional dark matter is therefore difScult to de- 
tect other than gravitationaly, despite its very consider- 
able flux on the Earth. 

Over the years, the idea that the dark matter might 
not be so weakly interacting, whether with itself or with 
ordinary matter, has been explored intermittently. In 
the late 1980s and early 1990s, the idea that dark mat- 
ter could carry ordinary charge was explored [H, [li S S] 
and severely constrained. Severe constraints were placed 
on the dark matter scattering cross section off ordinary 
matter With the exception of some windows at 

lower mass, the general conclusion was that the dark 
matter could not interact with ordinary matter except 
weakly (tr ^ 10~'^° cm^) (and much more weakly at typ- 
icaly WIMP masses of 10 — lO'^GeV). One clear excep- 
tion was that if the mass of the dark matter was suf- 
ficiently large, then the number density of dark matter 
particles, and hence their flux on any natural or artifi- 
cial "detector", would be too low to permit any useful 
constraints. Thus cTproton-DM < 10"^"^ cm2(m/Mpianck) 
is unconstrained by any known astrophysical or detector 
limits. (Interestingly, if m < 1 GeV, then most limits at 
high cross-section also fail for a variety of reasons.) In- 
teresting generic limits are available again only when the 
dark matter is sufficiently massive that its gravitational 
interactions start to affect galaxy dynamics. Spergel and 
Steinhardt revived the idea 0] of strongly self-interacting 
dark matter several years ago as a way to explain the ab- 
sence of central cusps in galaxy cores described by stan- 
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It has been suggested [l3| that the dark matter could 
be in the form of compact composite objects (CCOs here- 
after) - "strongly interacting composite macroscopically 
large objects which [are] made of well known (sic) light 
quarks (or/and anti-quarks)." Such objects, also known 
as QCD balls, are "formed from ordinary quarks [or an- 
tiquarks] during the QCD phase transition when (sic) 
axion domain walls undergo an unchecked collapse due 
to the surface tension which exists in the wall" [14|. An 
important prediction of [14| is that the minimum baryon 
number of a CCO for which its internal Fermi pressure 
renders it absolutely stable against the surface tension in 
the axion domain wall surrounding it during formation 
is Bqco — 10'^^. However, metastable CCOs may form 
with Baryon numbers as small as 10^*^, based on limits 
from the non-detection of neutral soliton-like objects by 
the Gyrlyanda experiments at Lake Baikal [3| ■ 

With the mass of the CCO Afcco ^ ^CCO GeV, tra- 
ditional limits on strongly interacting massive particles 
are clearly not applicable - the flux of CCOs is far 
too low to register in any detector, including most as- 
trophysical ones. For Bqco — 10^'^, CCOs impact the 
Earth a few times per year, the Sun perhaps 10** times 
per year, and a typical neutron star, just once every 10'°^ 
years. Individual impact events on Earth might be visible 
to cosmic ray detectors, but the instrumented area of the 
Earth is far too small. However, it has been highlighted 
that seismic shock waves resulting from the passage of a 
CCO through the Earth may already have been detected 

0- . . . 

Individual impacts of CCOs with the Sun are un- 
likely to be observable. Because the CCO bulk mat- 
ter consists of extremely dense superconducting anti- 
quarks/quarks, the majority of hadrons traversing the 
interstellar medium possess kinetic energies far below the 



2 



superconducting gap and are therefore unlikely to pen- 
etrate the ceo before being elastically scattered. The 
geometric cross-section of a CCO is given by 



<^cco — T^R, 
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where the CCO radius i?ccOj is determined by equating 
the Fermi pressure in the bulk matter with the pressure 
associated with the surface tension in the domain wall 
forming it. FoUwing the treatmncnt in 14], the typical 
radius is then given by 
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where c ~ 0.7 is related to the degeneracy associated 
with the massless degrees of freedom of each CCO and 
a ^ 10*^^^ GeV^ is the axion domain wall tension, which 
is constrained by axion search experiments. A CCO 
with Bcco = 10'^^ will therefore possess a typical radius 
Rcco ~ (10 — 100)/Ltm, and strike 
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nucleons on its way through the Sun, transferring ap- 
proximately 10^^~^^ GeV of energy, but only at a rate of 
1012-14 GeV cm-i , or lO^^-i^ergs-i. This is IQ-^i^-ie) 
of the solar luminosity. Nevertheless, this does suggest 
that if -Bcco < 10^^ then CCOs will be significantly 
slowed down, and therefore captured within the Sun. At 
Bcco — 10'^'^, a maximum of ~ 10^^ CCOs could have 
been captured by the Sun to date. It is difficult to see 
how to argue generically that the Sun has not captured 
this number of CCOs, since 10*** baryons represents less 
than IQ-i'^ of the Sun's baryon number. Furthermore, 
one expects the CCOs, being so heavy, if stable inside the 
solar environment, will settle to the centre of the sun. If 
the CCOs were point particles, they might form a black 
hole which would consume the sun [6| , but they are com- 
posite objects of much too low a density to form a black 
hole of such low mass. Instead, they are likely to either 
dissolve in the Sun, or combine into a single QCD ball 
of increasing size, until they reach the maximum stable 
mass, after which additional CCOs will indeed dissolve. 

Similar considerations apply to CCOs that strike neu- 
tron stars. Individual impacts while reasonably ener- 
getic, are not sufficiently so to be observed across the 
galaxy - even the complete annihilation oi a B = 10^"^ 
CCO releases only 10"^° ergs, typically over a (significant) 
fraction of a second. (The luminosity of the Sun, by com- 
parison, is 4 X lO'^^ ergs~i.) However, much of the energy 
release occurs deep within the neutron star, so that the 
most of the energy eventually goes into heating the neu- 
tron star. The rise in temperature of the neutron star as 
a result of such an event is less than 1 degree. 



One of the claimed attractions of CCO dark matter is 
that it could potentially explain several astronomical ob- 
servations which indirectly indicate the possible presence 
of dark matter. Chief among these are the 511 keV line 
signal from the centre of the galaxy and the apparent 
excess in background gamma-rays of energies approxi- 
mately in the range 1-20 MeV. 

Firstly, we consider the 511 keV gamma-ray line from 
the galactic centre measured using the SPI spectrometer 
of the INTEGRAL satellite This annihilation signal 
is thought to be produced by the process 
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Spectral analyses of the line signal strongly indicate 
that such annihilations occur following the decay of 
positronium atoms (25% as para-positronium and 75% as 
ortho-positronium) consisting of electrons and positrons 
traversing the interstellar medium (ISM) [l^ . It has been 
proposed that such atoms may form when low energy 
electrons traversing the ISM, interact with low energy 
positrons contained within the atmosphere of positrons 
(or positronsphcrc) predicted to surround each antimat- 
ter CCO @, HjIj since for CCOs in the favoured colour- 
flavour-locked (CFL) superconducting phase, leptons are 
prohibited from entering the bulk matter (see e.g. [S]). 
Each positronsphere is electrostatically bound to an an- 
timatter CCO by its extremely large net negative charge 
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resulting from a deficiency of strange antiquarks in a thin 
surface layer of the CCO bulk matter, owing to its fi- 
nite surface area [l^ . The positronium atoms arc then 
thought to survive until they decay, upon which they emit 
two 511 keV photons, which are claimed to contribute to 
the observed line signal. 

In [is'l, it is claimed that the observed rate of 511 keV 
photon production is "not in contradictions withobserva- 
tions for sufficiently large -Bcco"- This is obviously true 
since the annihiliation rate in the above scenario per unit 
volume per unit time, as a function of the distance, r, 
from the galactic centre is given by eq. (2) of Q as 
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where pB,nB are the energy and number densities of 
galactic baryons repectively, pdm is the energy density of 
dark matter particles (which in the present context are 
CCOs), and Uici. is the relative speed between colliding 
baryons and CCOs. Thus, we observe that a sufficiently 
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large value of Bcco will yield an annihilation rate con- 
sistent with the observations by INTEGRAL. 

However, the rate ^ assumes that every electron in- 
cident on a ceo will form positronium and ignores the 
relative suppression owing to the fact that such electrons 
will only be able to interact with the lower density re- 
gions of the positronsphere due to the increasing levels 
of electrostatic repulsion as it penetrates further into the 
positronsphere. In what follows we estimate the value of 
the suppression factor P associated with this effect. But 
firstly we must have a description of the electrical proper- 
ties of the positronsphere surrounding each CCO. To do 
so, we follow the treatment by Hu and Xu in [i^, which 
is largely influenced by the original calculations relating 
to strange stars by Alcock, Fahri and Olinto [25| . 

The standard calculation proceeds by assuming that 
the quarks and positrons near the surface of the CCO 
bulk matter are locally in thermal equilibrium; the rela- 
tionship between the charge density and electric potential 
is described by the classical Poisson equation. The ther- 
modynamic potentials fJj as functions of the chemical po- 
tentials fii (where i = u,d,s,e~^ for up, down and strange 
antiquarks and positrons respectively), the strange quark 
mass rris and the strong coupling constant ac can be 
found in the literature j25| . Chemical equilibrium be- 
tween the weak interactions involving the three antiquark 
flavours and positrons is maintained by the conditions 
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and their number densities are determined and related 
by the relations 
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where z is the altitude above the surface of the surface 
of the bulk matter. 

The quark number densities drop to zero for z > 0, but 
are not uniform for z < when one correctly accounts 
for the finite surface area of the CCO. This results in a 
significant depletion of s in a thin surface layer of the bulk 
matter (with a thickness of order 1 fm), resulting 

in a net negative charge of the bulk matter Qcco ) given 
by equation ^ 19]. The chemical potential /i can then 
be determined by ensuring that the surface pressure of 
the bulk matter is zero. 

In the widely-adopted Thomas-Fermi model, the 
positrons are approximated by a non-interacting Fermi 
gas, with the Fermi momentum of the positrons at 
each altitude being equal to the electric potential V. The 
positron number density is then given by 



V3 

ne+{z) = for z > 0. 
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Note that the thermodynamical potentials fli are de- 
fined for V — Q and for zero temperatures. This is cor- 
rected by substituting /i^ jii — QiV in ([5]), where q are 



the respective electric charges [2J], and by adopting the 
reasonable scenario where the temperature of the CCOs 
is much less than the potential energy at z = 13, 2^. 
Adopting such a scenario yields 
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where the complex term ns{V) can be derived from the 
chemical potentials fli, but we omit here for clarity. 
The boundary conditions for Eq. p^ are 
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where Vq is the electric potential in the deep core of the 
CCO, and for present purposes, has a value that is almost 
indiscernable from that of the surface potential V^- 

By integrating Eg. p^ in (— oo,V^] and [Vc-l-oo) re- 
spectively, and invoking the boundary conditions (|13p . 
corresponding expressions for the electric field E = 
—dV/dz are obtained. The surface potential Vc = V{z — 
0) is clearly a function of etc, but here we adopt the con- 
ventional value of Vc — 20 MeV, used widely throughout 
the relevant liteature (see e.g. [13, [111), and obtained 
by substituting the canonical CCO bulk matter density 
nq{0) = ne(0) ~ 9no into eq. ((TT|l . where uq — 0.15 fm^'^ 
is the characteristic nuclear density. 

Substituting Vc into the equations for E and then inte- 
grating them, expressions describing the electrical prop- 
erties of the CCO positronsphere for z > are ^ 



V{z) = 



V 



VcZ 



(14) 



^ We should note that the 1-D approximations ea. lll4l l- (ll6t are 
particularly appropriate for altitudes z < -Rcco i ^^'^ despite the 
fact that we utilise them in our calculations at larger altitudes 
it is unlikely that the deviations from an exact 3-D treatment 
would significantly alter our conclusions. 
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FIG. 1: Classical trajectory for an electron incident with the 
positronsphere of a CCO with an initial relative speed /3""', 
equal to 10~^ (top) and 10"'^ (bottom). 
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We can utilise equation (|15p to calculate the approxi- 
mate trajectory of an incident electron traveling along a 
radial vector using the classical equation of motion 



dt2 



eE(z) 
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Figure [T] displays the classical trajectory of an elec- 
tron approaching a CCO from 2 = 00 with initial speeds 
(3 = 10"^ and 10~^, which are the approximate ex- 
tremes of the range of speeds expected for the majority of 
elctrons traversing the ISM, given the rotational velocity 
of dark matter and baryons in the galaxy. Clearly, we ob- 
serve that such electrons only interact with positrons at 
altitudes z > 6 /.tm, which in fact, as mentioned above, is 
of the same order as the typical radius of the CCO itself. 

We also need to ensure that quantum tunneling effects 
do not significantly increase the rate of interactions be- 
tween positrons and electrons, owing to electrons pene- 



FIG. 2: Numerical solutions for the probability amplitude 
l-i/jp, of electrons approaching the surface a CCO with an 
initial relative speed /3"i'', equal to 10~^ (top) and 10~^ (bot- 
tom). The black (solid) vertical line corresponds to the classi- 
cal minimum altitude ^nJfn^', the magenta (dashed) and green 
(dot-dashed) vertical lines correspond to the altitudes z^"'" 
and z^^"^' respectively, defined such that 90% and 95% of the 
probability amplitude for < z < Zmfn^' lies at larger alti- 
tudes. The probability amplitude is normalised so that its 
sum for z < z^fn!* is equal to unity. 
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trating the classically-forbidden region of the positron- 
sphere associated with significantly larger positron den- 
sities. To do this we use the Numerov method to nu- 
merically integrate the non-relativistic time-independent 
Schrodinger equation from within the classically forbid- 
den region to larger altitudes. The results for the prob- 
ability amplitude ji/'P, normalised so that the summed 
amplitude for altitudes below the classical minimum al- 



titude 



class. 



is equal to unity, are displayed in figure [2l 



We clearly observe that the majority of the probability 
amplitude within the classically-forbidden region (i.e. al- 
titudes below the black (solid) vertical line) is extremely 
close to the classical minimum altitude, z^?^"- (95% and 
90% of which is contained above the magenta (dashed) 
and green (dot-dashed) vertical lines located at altitudes 
respectively). To demonstrate this quan- 



^95% 



and 



,90% 



titatively we calculate the following ratio 
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We find that R = 1.0104 for /3 = lO^^, and i? = 1.0022 
for /3 = 10~'^, indicating that the increased density ex- 
perienced by the tunnehng electrons, relative to that ex- 
perienced by purely classical electrons, is suppressed to 
within 1 percent of the density evaluated at the classical 
barrier z — zf^^f^^^-. Hence, we consider this to be ample 
justification for adopting the classical trajectories of elec- 
trons incident upon CCOs when calculating their rates 
of interaction with positrons residing in their positron- 
spheres. 

Because of the rapidly decreasing nature of the 
positron density with increasing altitude (as demon- 
strated by eq. ), we expect that the probability for 
forming positronium is significantly reduced relative to 
when the incident electrons are allowed to traverse un- 
hindered to the ceo surface. The value of this relative 
suppression P, will be a function of the CCO surface 
potential Vc, aswell as the initial speed and impact 
parameter b of the incident electron, relative to the CCO. 
The suppression factor P can be expressed as 

P = exp ^- Tle+{z[t]) Cre-e+^P,(Wrcl) Wrcl dt^ ■ 

(19) 

We can obtain an approximate lower limit on P by 
adopting a radial electron trajectory, which approxi- 
mately maximises the path integral of positrons, and 
by using a value of the parapositronium cross-section of 
(Te-e+^P^ (wioOwrci ~ TTfln /3™ c, which is a generous up- 
per estimate (see e.g. [21| for the relationship between 
CTe-e+^p, and /3™ ). 

Using the above method, we obtain P < 10^"^ for 
= 10-2 and P < IQ-^ for Pf"}- = 10~^. There- 
fore, even if we adopt the estimate for the 511 keV flux 
$511 from positronium formation proposed by 

$ ~ 10-3 cm-2 , (20) 

which is of the same order as that observed by the IN- 
TEGRAL satellite, if we now take into account the sup- 
pression on this flux dscribed above by multiplying it by 
P, we obtain 

/ 1033 \i/3 

<I>5ii = 10-^cm-2s-M-i^ P, (21) 

V ^GCO J 

hence. 



if 33 \ i / o 
10-s cm-2 ( for /3''"- = 10"^ 

\BccoJ 
/ 33 \ 1/3 
10-^ cm-2 f -g^j for = 10^3 , 

which is consistent with the INTEGRAL observations 
for Bcco < 10^^ for = 10^^ 2 ^ reconcil- 
able for electrons incident with speed = 10~3^ even 
when using the fore-mentioned experimental lower limit 
Bcco = 10^°. 

We should also note that in much of the literature 
relating to strange stars composed of similar colour- 
superconducting matter (e.g. [l^H^), it is proclaimed 
that such stars possess positron/electrospheres which 
obey equations identical to eqs. (fT4|) -(fT6 |) . but are only 
several 1000 fm in depth, compared to the infinitely ex- 
tending atmospheres adopted in the present study. If 
we were to adopt such truncated positronspheres, the 
charge associated with the bare antiquark matter would 
not be entirely shielded, yielding CCOs (plus positron- 
spheres) with a net electrical charge. Even if this charge 
is small enough to be consisent with observations relat- 
ing to charged dark matter, the CCOs will still possess 
electric fields which repel incident electrons as discussed 
above. Furthermore because of the reduced screening as- 
sociated with these truncated positronspheres these fields 
will be more intense at large, positron-deficient altitudes. 
Thus not only will the incident electrons be repelled at 
altitudes larger than those calculated above, based on 
the above calculations they will be repelled far above 
the truncated positronsphere, and the interaction rate 
with positrons will be significantly less than that calcu- 
lated above, yielding a significantly smaller value of the 
suppression factor P. Hence, the above scenario should 
be considered to be 'optimistic' as far as explaining the 
511 keV line signal with CCOs is concerned, and we will 
continue to adopt this scenario throughout. 

It has also been proposed that the excess of gamma- 
rays detected by COMPTEL at energies ~ 1 - 20 MeV 
[26| can be naturally explained by the photons pro- 
duced in the (non-resonant) direct annihilation process 
involving^interstellar electrons and positrons surrounding 
CCOs ^7|. The appeal of such an explanation is that the 
energy at which the excess occurs is similar to the Fermi 
energy of the positrons located at the CCO surface, and 
hence direct annihilations involving such positrons would 
yield photons which may potentially contribute to this 
excess. 



^ However we should like to point out that such speeds for the 
majority of electrons traversing the ISM are optimistically high 
in conventional models describing the rotational velocity profile 
of the galaxy. 
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However, in it is assumed that the incident elec- 
trons can penetrate the positronsphere "unhindered" , 
right down to the surface of the CCO, and directly anni- 
hilate with positrons whose Fermi energy (which at any 
altitude is equal to the electric potential energy V{z), 
see eg. p^ ) is approximately ~ 20 MeV. From the 
above discussion we see that the penetration of interstel- 
lar electrons to altitudes below the classical minimum 
altitude z^f^"-, which is of order 10^ fm for /3j^"J- = 10"^ 
and 10^ fm for — 10~^, is exponentially suppressed 
due to the rapidly increasing electric potential as it ap- 
proaches CCO surface. 

Consequently, the electrons will only be able to directly 
annihilate with positrons at 2; > z^"**' with Fermi ener- 
gies hf{z < zSlfn" ) < ^(^nlfif ) = where T^"^'- is 
the kinetic energy of the incident electron, and equal to 
~ 3 X 10"^ MeV for /3™- = IQ-^ and ~ 3 x 10"^ MeV 
for /3™- = 10^2, i.e. in both cases < 1 MeV. 

There is a finite probability that such electrons may 
quantum tunnel through the classical potential barrier 
to the surface of the CCO, and here we estimate an up- 
per limit for which as follows. The solution for the elec- 
tron probability amplitude |'(/'(z)^| locally around a given 



altitude z < z^ira!*' is of the form 



|?/'(z)p cx ex.p[K{z)z], k{z) 



2m (Viz) - T) 



1/2 



(22) 

Since V monotonically increases with decreasing z, the 
value of the exponent k also increases with decreasing z, 
so that we may say 



1^(0)1^ 



IV'(^'>^mfn")l 



< 



-exp[ 



-z'k{z')] , (23) 



where the LHS is the probability of an incident electron 
located at z^f^^- to quantum tunnel to the CCO surface, 
and the RHS is obtained by extrapolating the probability 
density at z' using a fixed exponent. Here for simplicity 
we use z' = z^^^" and obtain 



In 
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< 



-388 for = IQ-^ 

-1808 for = 10-3 



which effectively means that the probability for an elec- 
tron to tunnel to positrons which possess Fermi energies 
large enough for their direct annihilations to contribute 
to the MeV gamma-ray excess is vanishingly small. 

However, electrons resulting from pair production, due 
to the extremely strong electric fields near to the CCO 
surface (27I [28j , may directly annihilate with positrons 
in the surrounding positronsphere at sufficiently small 
altitudes in order to produce MeV gamma rays which 
may contribute to the MeV excess. To calculate such 
fluxes requires a knowledge of the thermal structure of 
each CCO which lies outside of the scope of this study 



We should also mention that positronium atoms 
which may form close to the CCO surface involving 
electrons either resulting from pair production, or those 
undergoing extremely improbable tunneling events, 
would likely be quickly reionized by the intense static 
electric field originating from the bulk anti-quark mat- 
ter. We expect that positronium atoms would quickly 
destabilize or have its initial formation suppressed 
when a static electric field exceeding strengths of order 
13.6 eV is present. Substituting Vc = 20 MeV 

into eq. (|15p . we observe that E < 13.6 eVA ~^ for 
z > 3.34 nm. Hence, we expect that Positronium 
formation will be suppressed inside this radius due to 
this effect, which for tunneling electrons will be in addi- 
tion to the extreme levels of suppression discussed above. 

Conclusions:- In this study we have considered the pos- 
sibility of compact composite object dark matter, com- 
posed of colour-superconducting di-antiquark pairs, as 
a solution to several unexplained astronomical observa- 
tions. Firstly, we considered the 511 keV line signal 
emerging from the galactic centre and observed by INTE- 
GRAL, and secondly the excess of background gamma- 
rays with energies approximately in the range 1-20 MeV, 
and observed by COMPTEL. 

It was proposed that the 511 keV excess could be ex- 
plained by the decay of parapositronium atoms formed 
from interstellar electrons and positrons electrostatically 
bound in a positronsphere surrounding each CCO. It was 
also proposed that the COMPTEL MeV excess could be 
produced from the direct annihilation of interstellar elec- 
trons and high energy positrons located near to the sur- 
face of the bulk matter of the CCO. 

In such proposals it was assumed that the incident 
electrons could travel unhindered to the CCO surface, 
however in this study we have deduced that insterstellar 
electrons approaching these CCOs with typical speeds 
will be unable to reach altitudes below approximately 
10^ fm (compared to zero). Consequently, the rates of 
positronium formation and the rate of direct annihila- 
tion at the CCO surface are significantly reduced; in the 
case of positronium formation the relative suppression 
P < 10~^, and the probability for direct annihilation at 
the CCO surface is exponentially suppressed to zero since 
such annihilations can only occur by incident electrons 
which tunnel through the potential barrier. 

Hence, in conclusion we find that the respective 
contributions of these processes to the fore-mentioned 
excesses are significantly reduced relative to previous 
estimates, to the extent where a CCO explanation of 
INTEGRAL and COMPTEL excesses is significantly 
less motivated. 
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